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Abstract 

Cascading gauge theories of Klebanov et. al. provide a model within a framework of 
gauge theory/string theory duality for a four dimensional non-conformal gauge theory 
with a spontaneously generated mass scale. Using the dual supergravity description 
we study sound wave propagation in strongly coupled cascading gauge theory plasma. 
We analytically compute the speed of sound and the bulk viscosity of cascading gauge 
theory plasma at a temperature much larger than the strong coupling scale of the the- 
ory. The sound wave dispersion relation is obtained from the hydrodynamic pole in the 
stress-energy tensor two-point correlation function. The speed of sound extracted from 
the pole of the correlation function agrees with its value computed in |hep-th/0506002 
using the equation of state. We find that the bulk viscosity of the hot cascading gauge 
theory plasma is non-zero at the leading order in the deviation from conformality. 
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1 Introduction 



The correspondence between gauge theories and string theory of Maldacena [1,2] has 
become a valuable tool in analyzing near-equilibrium dynamics of strongly coupled 
gauge theory plasma [3-10,13,11,12,14-17]. The research in this direction is primarily 
motivated by its potential application for the hydrodynamic description of the QCD 
quark-gluon plasma believed to be produced in heavy ion collision experiments at 



Previously, the dual supergravity computations [8] were shown to reproduce the 
expected dispersion relation for sound waves in strongly coupled Af = 4 supersymmetric 
Yang-Mills (SYM) theory plasma 



where v s , T], s, ( are the plasma sound speed, shear viscosity, entropy density, and bulk 
viscosity correspondingly. Conformal symmetry of the M = 4 gauge theory insures 
that 

v '=-j$> ^ = - ( 2 ) 

Non-conformal gauge theories, and QCD in particular, is expected to have nonvanishing 
bulk viscosity. Additionally, the presence of a scale in a gauge theory breaks conformal 
invariance responsible for speed of sound expression in Eq. (J2J). To see this note that in 
conformal theories the energy density and the pressure are related as follows e = 3P, 



The first computation for the speed of sound and its attenuation in strongly coupled 
non-conformal four dimensional gauge theories from the dual supergravity perspective 
were reported in [17]. Specifically, the model considered there was a mass deformation 
of the M = 4 Yang-Mills theory by giving identical masses mj to bosonic components 
of the two SYM chiral multiples and identical masses m/ to fermionic components of 
the same chiral multiples. Generically m/ ^ mj, and the supersymmetry is completely 
broken. For a special case m/ = m\, such a deformation leaves eight supersymmetries 
unbroken, and the model is usually referred to as M = 2* gauge theory [21-24]. Us- 
ing finite temperature gauge/gravity correspondence for this mass deformed M = 4 
SYM [25] the speed of sound and the ratio of the shear to bulk viscosity was found, 



RHIC [18-20]. 




(1) 



thus, vl = f = f . 
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respectively, 
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where ~ 0.9672, 01 « 8.001, and the ellipses denote higher order terms in rrif/T 
and rrib/T. From the dependence in Eqs. ©, it follows that at least in the high 
temperature regime the ratio of bulk viscosity to shear viscosity is proportional to the 
deviation of the speed of sound squared from its value in conformal theory, 

~ (5) 



V 

where k = 3tt(3^/2 « 4.558 for m& = 0, and k = 7r 2 /3^/16 4.935 for m/ = 0. 

In this paper we extend analysis of the transport properties in strongly coupled 
non-conformal gauge theory plasma to theories with spontaneously generated mass 
scale. More precisely, we analytically compute the speed of sound and its attenuation 
in 'cascading gauge theories' [26-28] (see [29] for a recent review) at temperature much 
higher than the deconfinement and the chiral symmetry breaking scales of the theory. 
The relevant finite temperature deformations of the theory were discussed in [30-32]; 
the holographic renormalization of the cascading gauge theories was explained in [33]. 
Moreover, in [33] the equation of state describing hot cascading gauge theory plasma 
was obtained from which the speed of sound was predicted to be 

BV — 1 9 
s de § 3 91n£ V ; 

where A is the strong coupling scale of the theory and ellipses denote subdominant 
terms for T ^> A. Here, we extract the dispersion relation ((TJ for the cascading gauge 
theory plasma from the pole of the thermal two-point function of certain components of 
the stress-energy tensor in the hydrodynamic approximation, i.e., in the regime where 
energy and momentum are small in comparison with the inverse thermal wavelength 
(uj/T <^ l,q/T 1), and at high temperature T ^> A. The latter computation is 
equivalent [16] to determining the dispersion for the lowest quasinormal mode in the 
corresponding black brane geometry (in our case [32]). The general prescription for 
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computing the quasinormal modes introduced in [16] has been applied recently to a 
variety of gauge/gravity duality examples [34,17,35]. 

As the computations are rather technical, we begin by summarizing our results 
in the next section. In section 3 we review five- dimensional effective supergravity 
description dual to chirally symmetric phase of the cascading gauge theory [33] and 
the explicit analytic construction of the background black brane geometry [32] dual to 
thermal cascading gauge theory at temperature much larger than its strong coupling 
scale. In section 4 we study fluctuations of the corresponding black brane geometry 
dual to a sound wave mode of the cascading gauge theory plasma. We introduce gauge 
invariant fluctuations and obtain their equations of motion. These equations of motion 
are valid beyond the hydrodynamic approximation, and for arbitrary temperature (as 
long as it is higher than the chiral symmetry breaking scale of the cascading gauge 
theory). In section 5 we derive and solve fluctuation equations in the hydrodynamic 
limit and at temperatures much higher than the cascading gauge theory strong coupling 
scale. Imposing Dirichlet condition on the gauge invariant fluctuations at the boundary 
of the background black brane geometry determines [16] the dispersion relation for the 
lowest quasinormal frequency ((H). Using the universality of the shear viscosity to 
entropy density ratio in strongly coupled gauge theory [11,12,15], computed dispersion 
relation can further be used to evaluate the ratio of bulk to shear viscosities. Some 
computational details are delegated to Appendices A, B and C. 

2 Summary of results 

Cascading gauge theory at a given high-energy scale resembles M = 1 supersymmetric 
SU (K*) x SU(K* + P) gauge theory with two bifundamental and two anti-fundamental 
chiral superfields and certain superpotential, which is quartic in superfields. Unlike 
ordinary quiver gauge theories, an 'effective rank' of cascading gauge theories depends 
on an energy scale at which the theory is probed [30, 32, 33] 

K* = K,(E) ~2P 2 In j, £>A. (7) 

At a given temperature T cascading gauge theory is probed at energy scale E ~ T, and 
as T > A, K*(T) 3> P 2 . In this regime the thermal properties of the theory [32,33] 
are very similar to those of the M — 1 SU (K*) x SU (if*) superconformal gauge theory 
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of Klebanov and Witten [39], with 

& cascade = (8) 

A* 

being the deformation parameter. Clearly, at the temperature increases, ^cascade 

be- 
comes smaller and smaller. From the equation of state for the cascading gauge theory 
one obtains the speed of sound in cascading gauge theory plasma as 



1 4 

V ° = 3 ~ 9 6cascade + °^cade) ■ ( 9 ) 

In this paper we compute the dispertion relation for the lowest quasinormal mode 
in the black brane geometry holographically dual to thermal cascading gauge theory. 
At high temperature, 5 cascade < 1, we find (see Eqs. (fTIjl. (fTHJ) . (j9~lj). dTTOJ) ) 

Ufa) = -J= - \ Scascade^j Q ~ « ^ ^ + jj<5 C aW e ^ + O (j^, ^ascad^j ' ( 10 ) 

from which we precisely reproduce the speed of sound extracted from the equation of 
state (jOJ). Furthermore, comparing Eqs. $1} and (fTTH) and using the universal result 
[11,12,15] 

?7 1 , s 

S 47T 

we arrive at the 'phenomenological relation' for cascading gauge theory plasma similar 
to © 

C o f .,2 1 \ , ml X2 



vi--)+OUl scad£ j. (12) 
Most importantly, it appears that phenomenological relation 

(13) 



£~-o(i)x ^ 




2 1 
f . - - 



is a robust prediction for hot strongly coupled non-conformal gauge theory plasma 
no matter whether scale invariance is broken explicitly (by masses as in Eq. (jSJl) or 
spontaneously (by a strong coupling scale as in Eq. (fT2|) ). As such, we expect it to be 
of relevance to real QCD quark-gluon plasma. (Note that the result f!13|) appears to 
disagree with the estimates ( ~ 77 {v 2 s — 1/3) 2 [36,37], later criticized in [38].) 
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3 Effective actions and equations of motion 

Effective 5d action describing supergravity dual to cascading gauge theories is given 
by [33] 

4/ d i V7 -t,V7M-t>\ D 2n „-* / V M i\ V i\ e 



Ss = tsgi y M5 vo1 ^ 1"^° - 2 V ** V "*J - «^(^^- + m 



1 A' 2 



2 fixfi^ 



where Rio is given by 

R w = R 5 - 2^ l g Xv [ VaV^i ) - 8^V( V A V,fi 2 



(14) 



- Ag Xu [2 il^ 1 ^ 1 V A ^i V^ 2 + 3 {if V A fi 2 V^ 2 J ( 15 ) 
+ 24 ^2 2 - 4 0| ^2 4 > 
with R 5 being the five dimensional Ricci scalar of the metric 

dsl=g^(y)dy»dy», (16) 
and G*, is the five dimensional effective gravitational constant 

G 5 = ^. (17) 

VOiyl,! 

We find it convenient to rewrite the action ()14j) in 5d Einstein frame. The latter is 
achieved with the following rescaling 

Qitu -> 9nv = fii /3 ^2 /3 9nu ■ (18) 

Further introducing 

fi x = e^ 4 ™ , fi 2 = e /+ ™ , (19) 
the five dimensional effective action becomes 

* = / voU {* 5 - f - 200*0' - - ±- 2{ 8Kfe-~ 



-V). 

(20) 



where we defined 



16, „ 16 , ln n t 28 j. . . I „ 40 , 

V = -2Ae-^ f ~ 2w + Ae~^ f - 12w + p 2 e ^ f+Aw + -K 2 e~- f 

2 



From Eq. ()20|) we obtain the following equations of motion 

3 



= D/ + ^e^- 4 -(^) 2 



3 dV 

80 a/' 



(21) 



(22) 



o = r> + JL e -*-v- to (M) 2 - - ^ 

1 40P 2 V ; 40 dw 



(23) 



= D*+4pJ e 



-$-4/-4u;^^2 



9P 
3$ 



(24) 



= □ A: - dKd(<$> + 4/ + 4iu) - 2P 2 e q,+if+4w 



dV 
~dK 



(25) 



(26) 

3.1 Black brane background geometry 

Taking the black brane metric ansatz 



ds^ = —c\ dt 2 + c 2 , dx 2 + c 2 dr 2 



(27) 



and assuming that all matter fields {/, w, K} depend on the radial coordinate r only, 
the background equations of motion are 



=/" + /' 



ln 



CiC| 

c 3 



, 3 -$-4/-4^/rW\2 _ _3_ 2&P 

+ 80P 2 1 j 80 3 <9f' 



:28) 



=«/' + w' 



In 



C 1 C 2 
C 3 



1 1 rW? 

+ 40P 2 1 j 40 3 dw 



(29) 
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=$" + $' 



111 



C 1 C 2 
C 3 



+ _l_ e -.-4/-4» ( ^ )2 _ ci 



<9P 
9$ 



(30) 



=K" + if' 



In 



ci4 

c 3 



- AT'[$ + 4/ + Aw}' - 2P 2 c 2 ! e* +4/+4 < 



dP 
dK 



(31) 



-ci' + c[ 



ln^ 

c 3 



+ 3^, 



(32) 



=c 2 + c 2 



111 



C 1 C 2 

c 3 



+ ^c 2 c 2 3 V 



Additionally, there is a first order constraint 



(33) 



on i 

= ($0 2 + y(/') 2 +40K) 2 + ^e-*- 4 /- 4 -(ir') 2 -12[ln c 2 ]'[ln Cl c 2 ]' - 2c 2 P . (34) 

Notice that Eqs. (PHj) - (jM|) are equivalent to the set of equations derived in [31,32] 
provided we identify 



ci = e 



c 2 = e 



c 3 = e 



g /-*+6» ^ / = (35) 



where {x, y, z} are correspondingly {x,y, z] of [31,32]. We would need solution to 
Eqs. (128)1 - ()34j) to leading order in P 2 . This was originally done in [32], but we repeat 
the main steps to set up conventions. We find convenient to use a new radial coordinate 



x = 



c 2 



(36) 



In terms of x Eqs. fl28)) - (}3~4")) become 



=c: 



" " (4) 2 - -A + 4y(/') 2 + f K) 2 + \W- — • 1 iJ - Ur < 



C2 



12P 2 



(37) 



=f" + -f 4 



x 80P 2 



(iT) 2 + 



<91nP f 9 



df [20 



^(M) a + ^M-?K) 2 



4 



160 v ; 2 W ; 320P 2 



e -*"4/-4«'(A")2| , 



(3f 
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(39) 

=*» + ^ + ^e—/--(^) 2 + ^ {l2([lnc 2 ]') 2 + Vc 2 ]' - 20(,/) 2 

- \{&? ~ f iff ~ ^- $ - 4/ - 4w (^) 2 } , 

=K" + pus - $ - 4/ - H'if + P 2 e* +4/+4 "{24([lnc 2 ]') 2 + — [hic 2 ]' 

- 4ok) 2 - ($r - f (/o 2 - ^p- 2 e-*- if - 4w m 2 } ■ 

Notice that Eqs. (|37|) - (j41j) have an exact scaling symmetry 

c 2 ^Ac 2 , /-/, w^w, K^K, (42) 

for a constant A. We will see later that this symmetry has an extension for the fluc- 
tuations as well. Physically, this symmetry corresponds to choosing a reference energy 
scale. 

To order 0(P 2 ) the solution to Eqs. (ESJ-flHJ) takes form [32] 



(43) 



w = ^(x), $=L-S(x), K = K* + P 2 k(x), 
where a is a constant nonextremality parameter, and 

e = ^(l-ln(l-x 2 )), K = -iln(l-x 2 ), 

K 7T 2 1 1 1 

C = -j^^ho^zon + — - 2 dil °g( x ) + 2 dil °s( 1 + x ) " 2 lnxln ^ ~ x ) ' ( 44 ) 
^ 3(1 + x 2 ) / M , N vr 2 \ , i 1 ^ 



1 = -8^3^ ^ilog(l - x 2 ) - - ) + -(1 - ln(l - x 2 )) 
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Furthermore, ip satisfies the linear differential equation 



°^" + ^'-(T^-!o(T^r < 45 > 

with the boundary condition 

i> = ^horizon + C(X 2 ) , ^ = " ^) + O ((1 - X 2 ) 3/2 ) , (46) 

where the second boundary condition will uniquely determine Vw^on- 

Finally, the (exact in P 2 ) temperature of the nonextremal solution is given by 



(2ttTY 



(47) 



4 Fluctuations 

Now we study fluctuations in the background geometry 

9 '/mu * 9 ' i-lu ")~ hfiv j 

w — > w + Sw , (48) 
$ -> $ + 5$ , 

x -> a: + 5^ , 

where {g^, f, w, $, K} are the black brane background configuration (satisfying Eqs. 
(P8 ]) -(P1 |) ). and {hfj, v , Sf, Sw, 5$, SK} are the fluctuations. We choose the gauge 

h tr = h Xir = h rr = . (49) 

Additionally, we assume that all the fluctuations depend only on (t, x%, r), i.e., we have 
an 0(2) rotational symmetry in the X1X2 plane. 

At a linearized level we find that the following sets of fluctuations decouple from 
each other 

{h XlX2 } , 

{h XlXl h X2X2 } , 

{h tXl , h XlX3 }, (50) 
{htti h aa = h XlX1 + h X2X2 , h tX3 , h X3X3 , Sf, Sw, S<&, SK} . 
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The last set of fluctuations is a holographic dual to the sound waves in cascading gauge 
theory plasma which is of interest here. Introduce 



h tt 


= c \ h tt = e - iu>t+iqX3 


4 


H tt , 


h tz 


=c 2 2 htz = e~ lujt+iqX3 


4 


Htz , 


foaa 


_„2 i _ -iujt+iqx 3 
— L 2 lb aa — c 


4 


H a a i 


' ^ZZ 


=r 2 h = p -iut+iqx3 
*2 zz 


cl 


ZZ } 


Sf 


_^-iu)t+iqx 3 jz 






Sw 


_^-iu)t+iqx 3 q 






5$ 


_ e -iujt+iqx 3 






5K 


g— iwt+iqx3 f£ 








haa h zz , Hit 




H a a H zz 



(51) 



where {H tt , H tz , H aa , H zz , J 7 , Q,p, JC} are functions of a radial coordinate only. Ex- 
panding at a linearized level Eqs. ([22 ft - (|26 p with Eq. ()48|) and Eq. (joT|) we find the 
following coupled system of ODE's 



o =K + K 



2 S "i ' 

In^M 

C 3 



H' u [ln Cl ]' 



r 2 / r 2 



<f ^ H tt + uf Ha + 2coq H t 



tz 



dP dP dV dV 



K 



(52) 



=H't' z + H' tz 



In 



Cic 3 



+ 4^1 H a 



(53) 



=H'L + ff' 



In 



cic| 
c 3 



C l V C 2/ 



4 2 

+ -cl 
3 3 



9/ 



(54) 



=H' + i/^. 



In 



c l c 2 
C 3 



(H' aa -H^ [lnc 2 ]' 



c 



+ I cj z + 2cug F t , + q 2 ^(H u - H, 



(&p &p &p &p 



dw 
11 



(55) 



=T" + T 1 
3 



In 



31 1 

c 3 



+ - 2 r [h. 



U V _l_ C 3 ( ' , ,2 „2 C 1 



1 



J 7 



2 /W W W 0^ 

80 ° 3 V 9/ 2 <9/<9u> 9/9$ P d/dtf 



(56) 



p2 



=ft" + ft' 

--4 

40 3 
1 



In 



ci4 



, + ^'[^-^] , + |(- 2 -g 2 |) n 



c 3 

/ d 2 V r d 2 V „ d 2 V 
T H Q H 

\ dwdf dw 2 dwd§ 



P 



d 2 V 
dwdK 



JC 



(57) 



=p" + p' 

if. 



In 



c l c 2 
C 3 



c 



r 2 

4 



p 



d 2 V r d 2 V n d 2 V d 2 V 

^+ P + 



3 \d$df " ' d<$>dw ' <9$ 2 r ' d$dK 

'^ 2 1 K'K! 



(58) 



-<E>— 4f—4w 



=/c" + /c' 



In 



C 1 C 2 
C 3 



2 p2 e $+4/ + 4«, c 2 

<9P 



z c x \ c 2 / 

/ d 2 P 

) 



( d 2 V d 2 V d 2 V d 2 V 



dKdw 



dKd<$> 



OK 2 



(59) 



+ ^(p + 4F + 4ft) - /C'[$ + 4/ + 4iu]' -K'\p + 4F + 4ft]' , 

<7i\ 



where all derivatives dV are evaluated on the background geometry. Additionally, there 
are three first order constraints associated with the (partially) fixed diffeomorphism 
invariance 



=u H', + 



In 



£2 



HiA+ql H' tz + 2 



In 



C2 
Cl 



if, 



tz 



(60) 
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=q K - 



In 



2l 

Cl 



H tt 



U H' t z ~ 1 Ha 



(61) 



=[lnc lC 2]'^ - [lnc*]' ^ + °\ (uo 2 H u + 2uoq H tz + q 2 °\ (H tt - H a 



c 



, ,'dP r &P ^ &P dV „ 
v <9/ <9u> (9$ <9if 



(62) 



We explicitly verified that Eqs. (fH2"j) - (f5HJ) are consistent with constraints (jOOJ)— ^H2J) • 
Introducing the gauge invariant fluctuations 

Zh =4— H tz + 2 H Z g — ff fl a 1 7,—. — + 2—^ H tt , 



a;^ c 2 C2 



Z f =T 



[hT^] 7 Haai 



w' 



H n 



[lnc 4 ]' aa ' 



(63) 



[lnc 4 !' aa ' 



J A" 



-2J 

K' 



H n 



[lnc 4 ]' aa ' 



we find from Eqs. (fo^j) - (jo^j) . (|5UJl - (J52)l . decoupled 1 set of equations of motion for Z's 



—AhZ" h + BhZ' h + CrZh + DnZf + EhZ w + FhZ<$> + GrZk 



(64) 



=A/ZJ + fi/Z} + C,^ + DjZji + E f Z f + F f Z w + + if^ + J/Z* , (65) 

1 To achieve the decoupling one has to use the background equations of motion H28|) - (|34[l . i.e., the 
decoupling occurs only on-shcll. 
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(66) 




(67) 




w 



(68) 



where we collected connection coefficients {A..., • • • , K K } in Appendix |XJ Notice that 
Eqs. (JHIJ)-(|HEI) have an exact scaling symmetry (J4*2"|) provided the latter is supplemented 
with 



while keeping {Z H , Zf,Z w ,Z^>, Z K } invariant. 

5 Hydrodynamic limit, boundary conditions and small P 2 ex- 



We study now physical fluctuation equations ()64)) - (j68|) in the hydrodynamics approxi- 
mation, u — > 0, q — > with - kept constant, and to leading order in P 2 . Similar to 
the N = 2* computations [17], we would need only leading and next-to-leading (in q) 
solution of (}6~4*j) - (p)8*j) . The computations are greatly simplified with judicial choice of 
the radial coordinate. Choosing the radial coordinate as in ()36|). we find that at the 
horizon, x — > + , Z H oc x ±iuj ^ 27tT \ and similarly for Zf, Z w , Z®, Z K . The temperature 
T is given by Eq. (}4Tj) . Incoming boundary conditions on all physical modes implies 
that 



A LO 



q^ Xq, 



(69) 



pansion 



Z H (x) = x m z H (x) , 
Z<$>(x) = x~ irv z<s>(x) , 



Z f (x)=x-*°z f (x) t 
Z K (x) = x~ tw z K {x) , 



(70) 



where {zh, Zf, z w , z$, Zk} are regular at the horizon; we further introduced 



ID = 



2vrT 



q = 



2nT 



(71) 
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There is a single integration constant for these physical modes, namely, the overall 
scale. Without the loss of generality the latter can be fixed as 



z H {x) 



(72) 



In this case, the pole dispersion relation is simply determined as [16] 



z H {x 



(73) 



The other boundary conditions (besides regularity at the horizon and (J73J)) are [16] 



z f (x 
Let's introduce 



0, z w {x) 



, z<$> (x) 



0, z K (x) 



as-»l- 



(74) 



zh = [z u H>0 + P 2 4,o J+M (^4,1 + ^ 4,i ) > 

*/ = (4o + p 2 4o) + < q (4i + p2 4i) 



4,o + ^ 2 4,o J + < q (4,i + ^ 2 4,i J , (75) 
^$ — ^4,0 

+ P 2 zl ^j + i q (zl A + P 2 4 ^ , 

z K = ^4 >0 + p 2 4 i0 ^ + z q ^4 4 + P 2 4 4 j , 

where the lower index refers to either the leading, oc q°, or to the next-to-leading, oc q 1 , 
order in the hydrodynamic approximation, and the upper index keeps track of the P 2 
deformation parameter. Additionally, as we are interested in the hydrodynamic pole 
dispersion relation in the stress-energy correlation functions, we find it convenient to 
parameterize 



tv = ^=(l+(3 v P 2 ) -^l+/3 r P 2 J, (76) 

where the P 2 = coefficients are those of the M = 4 plasma, computed in [8], and /3 V , 
fir are constants which are to be determined from the pole dispersion relation (|73j) 



'H,0 



7 2 



(77) 



15 



Using the high-temperature non-extremal cascading gauge theory flow background 
parameterizations (J75J) . (|7Sjl. we obtain from Eqs. (Jo^ - fRjHj) 2 four sets of ODE's 
describing leading and next-to-leading order in the hydrodynamic approximation and 
O(P ) and 0(P 2 ) order in the deformation parameter. 

In the remaining part of this section each set is discussed in details. 

5.1 Equations and solution in 0(q°) and O(P ) order 
To order O(q ) and O(P ) we find the following set of equations 

0-\z° Y'l 1-3x2 L~° + | 32(2(23-^ + 1 + ^) Q 

U " [ ^ oJ + x(l + x 2 ) [ ^ oJ + 1 + a; 2 + (1 - x*)K* Zr >° ' 1 ^ 



-[4ol" + - [4ol' " (1 _ 8 X 2 )2 4,o + ^ [4,o]' + (1 _ X 2 )2 ^ 4,o , (79) 



o — [4,o] ' + ~ [4,o] ~~ _ 4,o + [4,o]' > (so) 



4,o]" + ~ [4,o] + IT - [4,o]' j 



0=[4,of + ^[4,or- (82) 

Notice that at this order, the fluctuations couple only through z° KQ) which by itself 
decouples. The most general solution of Eq. (J%2~|) is 



4,o = Ci+C 2 In x. (83) 
Regularity at the horizon and Eq. (|74j) imply that 

4,o = • (84) 

Vanishing of z° KQ = decouples all the remaining fluctuations. The pattern that we 
observe here extends to order 0((\ l ) in the hydrodynamic approximation and to order 
2 After rewriting them in radial coordinate 
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0(P 2 ) in the deformation parameter: 

■ the corresponding gauge invariant fluctuations couple only through z'^ 

■ z'k ... by itself decouples; 

■ subject to boundary conditions, the unique solution is z'^ ... = 0, resulting in decou- 
pling of {<% ..., zj , z'^ .„, fluctuations. 

Given Eq. (j84j). the remaining equations can be solved analytically. With boundary 
conditions (|72j) - (|7H) unique solutions are 

z h,o = 1 — x 2 , Zj = z° j0 = z^ q = . (85) 

For z° H q we reproduce the sound wave quasinormal mode in the near extremal D3-brane 
geometry [8, 17] to leading order in the hydrodynamic approximation. 

5.2 Equations and solution in O(q ) and 0(P 2 ) order 

Using Eqs. (f5i|). (JSHj) . to order C(q°) and 0(P 2 ) we find the following set of equations 



n r 2 „, , 1 ~ 3x 2 2 4 32(2(x 3 -x)^ + l + x 2 ) ^ 



" x(l + x 2 ) L ^ oJ " 1 + X 2 " (1 - x 4 )i^ ^° 

: 

3x 2 (l+x 2 ) V ^* 

r 2 l// r 2 i / 2 3/ij r 2 1/ ^ 2 

=[ ^° J + x [Zffil ~ (i-x 2 ) 2 ^° + tok: [ ^' oJ + (i - x 2 ) 2 k, 

+ —T^ ^TTT ( 40(x 2 - 1)V + 40x(4r? - k) - 7a; J , 



(86) 



(87) 



60(x 2 - l)x 3 K* 

r 2 i // r 2 i / 2 ^ r 2 1/ 

= K>,oJ + - i Z w,o\ - 7] Z w,0 + TJ7~ i Z K,o\ 

x {lx) bK. (g8) 

+ 30(, 2 -l)^ (20(x 2 -l)^ + 30^ + x), 

=[4,o]" + ^ [4,o]' + % [4,o]' + 3( ^ 2 _\ )g3 ^ ( 2 (* 2 " 1)C + *) - (89) 



=[4,of + I [4,o]' + 3(x2 ! 1)a , 3 ( ^ - ^ + x ) " ( 9 °) 
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Notice that Eqs. ([8lj J) -(pU |) are equivalent to Eqs. (|75 jl -(|82 )l of the previous section 
apart from 0(P 2 ) sources describing the deformation of the nonextremal cascading 
gauge theory geometry away from the near extremal D3-brane background (|43j) . The 
latter is precisely the reason that the pattern of coupling of fluctuations is the same as 
for P 2 = 0. 

Using Eq. , the most general solution of Eq. (J9*U)) is 

4,o = C i+ C 2 Ins. (91) 
Regularity at the horizon and Eq. ()74j) imply that 

4,o = • (92) 
Given Eq. (j92|) . the most general solution to Eq. (}86|) takes form 

4,o = Ci (V - 1) In a; - + C 2 (1 - x 2 ) - ^4 + 6K*p v J . (93) 

Regularity at the horizon implies that C\ = 0, and the boundary condition (f7H|) deter- 
mines 

ft = -3^, (94) 

in agreement with [33]. Finally, since z° H0 already satisfies Eq. (J72J), we must also set 
C 2 in Eq. (JDIIJ) to zero. Thus, we have 

4,o = • (95) 
Remaining fluctuation equations are discussed in Appendix FBI 

5.3 Equations and solution in 0{(\ l ) and O(P ) order 

Using Eqs. (jBU), (|%5*j). to order O^q 1 ) and O(P ) we find the following set of equations 



n _r y, , 1 3x2 r 1/ , 4 32(2(x 3 x)/t' + 1 + X 2 ) q , , 

° -^> lJ + x(l + x 2 ) [ ^< lJ + 1 + x 2 Zh > 1 + (1 - x 4 )iv* ^ ' (%) 



o =[4il" + \ m - tA^ 4i + 7^ [4J' + „ Lr, 4,1 , or) 
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=[z, 



w % l) 



X 



(1 



X 



2\2 



Z° 



(9* 



X 



(99) 



Notice that two sets of equations Eqs. (|78| ) -(J82 |) and Eqs. ()96| ) - ()100|) are equivalent. 
With boundary conditions (|72|l - (j73jl unique solutions are 

z h,i = z °f,i = z w,i = z %,i = z k,i = • (101) 

Again, for z° H1 we reproduce the sound wave quasinormal mode in the near extremal 
D3-brane geometry [8, 17] to order q 1 in the hydrodynamic approximation. 

5.4 Equations and solution in 0{(\ l ) and 0(P 2 ) order 

Using Eqs. (JHU), jS5), (1111)-®, j22), (EH, (HHH), to order 0(q x ) and C(P 2 ) we 

find the following set of equations 



r 2 ]H 1 ^X r 2 1/ ^ 2 i 32(2(x^ x)/t -|- 1 -f- X ) 2 

u = %,i 71 ; 2Y I z h,i\ + "r~ 2 ^ T\ 1^1? z k,i 

x(l + x z ) ' 1 + x z (1 — X 4 )-/!* 

8\/3 

~ 9-ftr, (x 2 + 1) x 2 

+ 6x (l - x 2 ) 2 (In a; + 3) f + x 2 {3K*p r - 6 - 21nx)Y 



(V) 2 (1 - x 2 ) 2 ((x 2 + 1) lnx + 3x 2 + 1) 



(102) 



r 2 i// r 2 1/ ^ 2 r 2 i/ ^ 2 

=l*/,iJ + - L*/,il ~ (1 _ X 2)2 + 77^ L**,il + (1_ X 2)2^ ^,1 

- m , 2 Jh^l? a ( 40 ^* ( X + i 1 ~ t4o]' (103) 
60 (ar — 1) (x 2 + 1) A*x d \ 



(l - x 4 ) (40r/(x 2 - 1) - 40xk - 7x + 160xr/)^ 
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10, 1J l^w.l 
X 



3QK*x 3 (x* - 1) 

(x 2 + 1) (20(x 2 - + 30xip + x) 



(1-z 2 ) 2 ^ 5K* L *^ 
(Wi^l -* 4 )[4,o]' (104) 



(x 2 



+ 1) (2(x 2 -l)C' + *)) 



(105) 



=[4,]" + \ [4,i]' - 9g3( y_ 1} ((* 2 - + *) ((* 2 + 1) Inx + 3x 2 + 1) . (106) 
Given Eq. (JBJ), the most general solution to Eq. ()1U6|) is 

4,i = Ci+C 2 lnar. (107) 
Regularity at the horizon and Eq. (J73)l imply that 

4,i = 0. (108) 
With Eq. (jl08J) . the most general solution to Eq. ()102|) takes the form 

4,i = Ci (V - 1) lnx - 2) + C 2 (1 - x 2 ) + ^ (W*/? r - 2) . (109) 
Regularity at the horizon implies that C\ = 0, and boundary condition ()73|) determines 

* = m- (110) 

Finally, since z° H0 already satisfies Eq. (|72^). we must also set C 2 in Eq. ()109|) to zero. 
Thus, we have 

4,1 = 0. (in) 

Remaining fluctuation equations are discussed in Appendix O 
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Appendix 

A Connection coefficients for Eqs. (EH)-(EE| 

A H =3u 2 c 2 c 2 2 c 3 cj ( -c 2 c 1 q 2 d 1 - 2c' 2 cjq 2 + 3tu 2 d 2 c 2 2 ) . (112) 



B H =uo 2 c 2 Ci y27uj 2 c 2 2 cic' 2 3 c 3 - 42c?gV 2 3 c 3 + 9c 2 2 Ciq 2 d 2 d 2 c 3 - 2c 2 2 c\q 2 'd 2 c\V 

, r> 3 2 2 / 3t> i o 2 2 2 / / / , a 32/2/ o 22/2/ n23 /2 / 

+ 2c 2 c x q c x c z P + 3c 2 c 1 g c 2 c x c z + <ociC x q c 2 c 3 — 3c2C x q c 2 c x C3 — 9u c 2 C\C 2 c 3 
+ 9u 2 clc 2 c 3 c'^j . 

(113) 



C H =c 3 cj 2 ^-3cj 2 c^c 2 c^ 2 CiC / 1 + 9u 4 c 2 3 c' 2 4 ~ 2clq 2 c{d 2 c 2 2 V + 3Qq 2 cfcf c[c 2 

- 2Aq 2 c\4 + \c\q 2 c\d 2 d x c\V - 2c\q 2 c\d 2 d\V + 3c\d 2 c\q^d x c 2 + Qc 2 3 d 2 c 4 x q A (114) 

1 o 2 /3 3/ ir 22 /2 2221 

— 12g C\C X c 2 c 2 — 15a; c 3 c 2 c x q c 2 I . 



D H =16c 3 g 2 c 2 i—C\d 2 + d x c 2 J ( — ^-c 2 (—Cic\Vu 2 c\ — 6d 2 u) 2 d x C2 + 6ci(d 2 ) 2 uj 2 

1 9V 3 N (U5) 

+ c\q 2 c\V)f - - 4 c l cl — {2d 2 c\q 2 - 3d 2 c 2 2 u 2 + c.c^q 2 )) . 

Eh =16c3g 2 c 2 I — c\d 2 + c x C2 J I — 10c2(—ciC 3 1 Vuj 2 c 2 . — Qd 2 uj 2 d x C2 + Qc\d 2 uj 2 

i op \ (116) 

+ c\q 2 clV)w' - -clc l — {2d 2 c 2 l q 2 - 3d 2 c 2 2 u 2 + c^q 2 ) . 
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F H = 16c 3 g 2 c 2 ^— cic' 2 + c^c^j ^— -c 2 (— cic 2 'Pcc> 2 c 2 , — 6c' 2 u 2 'c[c2 + 6cic 2 u 2 
+ c\q 2 clV)& - ± Cl cj^(2c' 2 clq 2 - 3c' 2 c 2 2 u 2 + c^g 2 ) j . 



+ 6 Cl c 2 V + c\q 2 clV)K' - i c 2 Cl |^(2 C2C 2g 2 - 3c 2 c 2 cu 2 + c^q 2 )) . 



A f =12c 2 c 3 c 2 c 2 -2c' 2 c 2 g 2 + 3c 2 c 2 w 2 - c x c 2 c x q 2 ) . (119) 



(117) 



(118) 



Bf =Yld 2 c 2 cx (^—2c l 2 c\q 2 + 3c' 2 c 2 cj 2 — Cic^g 2 ^ ^3c 3 Cic 2 — ojc'gCx + c 2 c[c^j . (120) 
C> = - clclclu 2 (JL^W + 2ca /'7>) . (121) 
£>/ = - c^cV (c lC ' 2 - c 2 c^ (^4^ + W^) . (122) 



E/ = ^ C 3 C 2 C l C 2 ( 2C 2C 2 <? 2 - 3C 2 C 2 ^ 2 + C^c'^ ^ + ^/'c? ^-6c 2 C 2 ^ 2 

+ 5c' 2 c 2 g 2 + c lC2 c' ig 2 ) ^ + 3 -fclcjc 2 J> 2 {-c 2 2 u 2 + c 2 g 2 ) 7> 



,3 , / 3 /2 _ < j,_4f_4 w) 22 „„ 2 2 2,on 222\/o'22 o ' 



2 2 
2 UJ 



+ c x c 2 (! x q 2 



(123) 
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Ff = ^c\c 2 2 c\c' 2 ^2c' 2 c 2 iq 2 - 3c' 2 c 2 2 u 2 + c.c^q^j + I8c 3 2 clcjc' 2 w' (-c\u 2 

+ cfq 2 ) ^ + Aclclclf [2c' 2 clq 2 - 3c' 2 c 2 2 uj 2 + c^q 2 ) ^ 

/ X 9 / (124) 

+ mcl4clfw'(-cW + c 2 iq 2 jV + ^ 2 clc,cV 2 K> 2 e-*-^ [2c' 2 c\ q 2 

— 3c' 2 c\uj 2 + CiC 2 c[q 2 



9 2 „2 „3 „/ „' „2 „2 q„/„2, ,2 , „ „ „/ „2^\ ^ 2 ^ ? , 9 : > ;> . /K / ,f 2 ,2 



Gf =7^c 1 c 2 c 3 c 2 \^2c 2 c 1 q - 3c 2 c 2 u + c x c 2 c x q j + — c 2 c 3 c lC2 $ \^c 2 u 

+ c\q 2 \ ^ + Aclclclf Ud 2 c\q 2 - 3c' 2 c 2 2 u 2 + c x c 2 c' x q 2 \ ^ 

/ \ 9 / (125) 

+ ^c 2 /'*' (-c^ 2 + c 2 , 2 ) P + ^c^c^'V*- 4 ^ [2c' 2 c\q 2 

o / 2 2 , / 2 

— dc 2 C 2 CJ + Ci 020^ 

^/ = - Y^^c^cac^'e-*- 4 ^ 4 - ^2c 2 c 2 g 2 - 3c 2 c 2 c 2 + c.c^q 2 ^ . (126) 



^ ,3„2„2„/ Z' o„'„2„2 g„/„2, ,2 , „ „ 



J / = ™ c 3 c i c 2 c 2 Sc^qg - 3c 2 c 2 u; + c x c 2 c x q 



20 3 1 2 2 y 2 11 2 2 

+ ^ c 2^c?c 2 K'e-*- 4 /- 4 - (-c 2 c 2 + c 2 g 2 ) ^ + ^c 2 /' ( 2 c 2 c 2 g 2 (127) 
- 3c 2 c 2 c 2 + Cl c 2 c' ig 2 ) H + ^ 2 c\c\c\fK<c-*-^ {-du 2 + c 2 g 2 ) V . 

A w = Yld 2 c z c\c\ ( -2c' 2 c\q 2 + 3c' 2 c 2 cj 2 - Cic^g 2 ) . (128) 



5 W =12c' 2 c 2 ci {^-2d 2 c\q 2 + 3c' 2 c 2 u; 2 - cic 2 c' 1 g 2 ^ ^3c 3 cic 2 - c 2 c' 3 ci + 02^03^ . (129) 
C w = - c\c\c\u? [2c 2 w'V + ^c 2 f^) • (130) 
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D w = - 4d4uj 2 ^cic' 2 - c 2 d^j ^2c 2 w'V + ^4^^ 



^^2^4 ^2c 2 c 2 g 2 - 3d 2 c 2 2 u 2 + c^q^j + 4c^cX( 2c' 2 c^ 

' „2, ,2 | „ „ / „2| ^ . Q J 53 ) <■/ / 2 ,2 , 2 2 \ 



- 3c' 2C2 V + desci^J — + 80^4/' (^- C2 V + c^J — 
+ ^c^/V (-c 2 c 2 + c 2 g 2 ) V + J^c^'V*" 4 ^ ( 2 c 2 c 2 g 2 



o / 2 2 , / 2 

3c 2 c 2 cj + Cic^g 



^ 3 / 2 2 /" o I 2 2 o / 2 2 i / 2 A ^ ^ , ' 

=— CgCjqCal 2c 2 c 1 g - 3c 2 c 2 w + c&ctf I ^ + 4c 3 c 2 c 1 w ( -6c 2 c 2 o; 

../ ..2.. 2 i / ..2 \ i I r;n, : '>, J, 2 ,,/2 / ,2 ,2 2, ,2 



+ 5c 2 Ci<f + c lC2 cW) — + 160c^c 2 c>' 2 ^-c 2 u; 2 + c*<f J P 



+ cic 2 qg 2 



y^c 2 c 2 c 2 (^2c 2 c 2 g 2 - 3c 2 c 2 ^ 2 + c^q^j + y^c^c^' ( -c 2 ^ 2 

4c 2 c 3 c 3 w ,f 2c / 2 2 _ 3c / 2 2 , c -jJ\ 



+ cfc-J ^ + Aci44w' [2d 2 c{q z - 3d 2 c z 2 u z + c x c 2 d x q j ^ 
+ 4cl4clw'& (-c 2 c 2 + c\q 2 ^ V + -^cW\d 2 K' 2 e-^-^ (V 2 c 2 g 2 



3c 2 d\uj 2 + C\C 2 d x q 2 



H w = - ^ 2 c\c\ C?> d 2 K'e-^-^{2d 2 c\q 2 - 3d 2 c 2 2 u 2 + c.c^q 2 ^ . 
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3 / \ d 2 V 

t 1 rt >'2tlc\!'2'2 o / 2 2 / 2 \ 

=^Q C 2 c 3 c i c 2 1 2c 2 Ci? -3c 2 c 2 u; + CiC 2 qg 1 dwdR 



+ ^ 2 4 c lclc' 2 K>e-*-^ (-c 2 c 2 + c 2 g 2 ) ^ + 4^«/ (2d 2 cW (136) 

- 3C 2 C^ 2 + C lC2C ' iq ^ H + A c 4 c 3 c 2 w /^ e -*-4/-4. ^ _ ^2 + ^ p 

A $ =12c 2 c 3 c 2 c^ ( -2c 2 c 2 g 2 + 3c' 2 c^ 2 - c^c'^ 2 ) . (137) 



=12c 2 c 2 C! ^— 2c' 2 c 2 g 2 + 3c 2 c 2 u; 2 — Cic^g 2 ^ ^3c3C!C 2 — c 2 c 3 Ci + 02^3^ . (138) 
(7$ = — c\c\c\uj 2 (6c' 2 ^ + 2c 2 &V^J . (139) 
D„ = - c\c lC tuj 2 (c lC ' 2 - c 2 c'^j (6c 2 ^ + 2c 2 &V^J . (140) 



F,, 12c 2 c 2 c^ 2 ( 2c 2 c 2 g 2 - 3c 2 c 2 cu 2 + Cl c 2 c' iq 2 j + 4^c^c 2 <l> / \2c' 2 c\q 2 

' 2 .2 J S\^ > 1 r>rtf\J>J>„1J f>( 2 ,2 1 ^2„2\ 



\ <9 2 P 



- 3c' 2C2 V + c lC2C ' ig 2 J — + 320c^c^/'(-c^/ + cfg 2 J ^ 
+ ^c^c 2 /'*' (-c 2 c 2 + c 2 g 2 ) 7> + ^c^'V*" 4 ^ ( 2 c 2 c 2 g 



3c 2 c 2 tJ 2 + Cic^g 2 I . 



F$ =12 C 3 C2 c 2 C 2 ( 2c' 2 c 2 g 2 - 3c 2 c 2 cu 2 + c^q^j -^-^ + Aclc\c\<$>' (2c' 2 cjq 2 



- 3c' 2 c z 2 uj z + c lC2 d iq z j — + 480444^' y-c z 2 u z + ciq z j ^ 
+ 160cjc 4 2 clw'& (-cjoo 2 + cjq 2 ^ V + ^ 2 c 2 2 Czc\c' 2 K l2 e-^- iw ^cjq 



0/22, / 2 
3c 2 c 2 cj + C\C 2 c x q 



(141) 



(142) 
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G $ =\2c\d 2 c\c\ [2d 2 c\q 2 - 3c' 2 c 2 2 u 2 + c^q^j |^ + 4^^$'c? ( -6c 2 cjuj 2 
+ 5c 2 c 2 g 2 + c lC2 c' ig 2 ) ^ + 4c^4c 2 $' 2 f-c 2 c 2 + c\q 2 \V 



+ 3c 3 c' 2 (-l^' 2 e-*- 4 /- 4 -c 2 c 2 - 4c 2 c 2 c 2 + 4c 2 c 2 g 2 ) ( 2 c 2 c 2 g 2 - 3c 2 c: 



(143) 



2 , , 2 

2 W 



+ CiCscig 2 



if* = - ■^(? 1 (%c 3 (Z 2 K'e-*- Af -* w [ 2c' 2 clq 2 - 3c' 2 c 2 2 uj 2 + c^q 2 ) . (144) 



/A d 2 V 



U = l2c 2 2 c\c 2 1 c l 2 ^2c' 2 c 2 l q 2 - 3c' 2 c 2 2 u 2 + c^c^q 2 ^ d ^ dR 

+ §- 2 44clc' 2 K'e-^-^ {-c 2 2 u 2 + c 2 g 2 ) ^ + Aclclc\V ( 2 c 2 c 2 g 2 (145) 
- 3c 2 c 2 c 2 + c lC2 cW 2 ) || + A^^'if'e-^ ^2 + c 2 q 2^ v 

A K = 12c 2 c 3 c 2 c 2 I -2c' 2 c 2 g 2 + 3c' 2 c 2 cj 2 - c x c 2 d x q 2 j . (146) 



B K =12c' 2 c 2 ci ^-2c' 2 c 2 g 2 + 3c 2 c 2 u; 2 - Cic^g 2 ^ ^3c 3 cic' 2 - c 2 c' 3 Ci + ^0^3^ 
+ 12c 3 c 2 c 2 c 2 ( Aw' + $' + 4/' j ( 2c' 2 c 2 g 2 - 3c' 2 c 2 cj 2 + Cic^g 2 ) . 



(147) 



= - c\c 2 



Iclu 2 ^12c' 2 P 2 ^e^ +4w + 2c 2 K'V^ . (148) 



D K = - cl Cl c A 2 u 2 (c lC ' 2 - c 2 c'^ {l2d 2 P 2 ^e" + ^ +Aw + 2c 2 K'V^ . (149) 



^ =48c 3 c / 2 i ; s: / c 2 c^ 2c' 2 c 2 g 2 - 3c' 2 c 2 2 uj 2 + c x c 2 d x q 2 . (150) 
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; iclc\c\K'[ 2c 2 c\q 2 - 3c 2 c 2 cj 2 + c^c^q 2 ) — 



+ 32clc 2 1 c 2 2 c' 2 P 2 e' s ' +4f+4w ( 20clc 2 q 2 f - 20c^ 2 /' - 9c' 2 c 2 2 u 2 + 3 Cl c 2 c' iq 



2 



+ 64c 2 g 2 ) H + ^c 2 c 4 X'/' ( -cW + c 2 g 2 ) 7> . 



=A8c3c' 2 K'clc 2 l ( 2c 2 c 2 g 2 — 3c 2 c 2 a; 2 + c^qg 2 



H K =24djc 2 c 2 c 2 P 2 e* +4/+4w ( 2c 2 c 2 g 2 - 3c 2 cjuj 2 + Cl c 2 c' iq 2 1 ' 9 ' 



dwdK 

— 

+ 96djc 2 c 2 c 2 P 2 e* +4/+4w f 10c 2 c 2 g V + 2c 2 c 2 g 2 - 10c^ V - 3c' 2 c 



2 2 
2 u; 



+ c.c^q^j ^ + mOclciclK'w 1 [ -c\uj 2 + c\q 2 )V . 



Ik —12c 3 c 2 K' c 2 c 2 ( 2c' 2 c 2 g 2 — 3c 2 c 2 a; 2 + Cic 2 c[q 2 



J K =24clc 2 2 cjc' 2 P 2 e' s ' +4f+4w ( 2d 2 c\q 2 - 3c' 2 c 2 2 u 2 + Cl c 2 c[q 



_8 2 V_ 

d$dK 



dV 

+ Aclclc\K' ( 2c 2 c 2 g 2 - 3c' 2 c 2 w 2 + c x c 2 d x q 2 ) — 



+ 24c^ 2 c 2 c 2 P 2 e ,I>+4/+4w ic lC2 c' iq 2 - c^V + c 2 c 2 gV + 2c 2 c 2 g 2 
- 3c 2 c 2 c 2 ) H + AclclclK'V [~c\u 2 + c 2 g 2 ) V . 
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K K =24clc 2 2 c 2 1 c' 2 P 2 e* +if+ ' lw \2c' 2 c\q 2 - 3c' 2 c 2 2 u 2 + c x c 2 d x q 2 

l/ i33r^'2(p/22,r-/22| / 2 1 

+ 4C3C2A q I -6c 2 c 2 u + hc 2 d{q + Ci 02(^3 J 



d 2 V 
OK 2 



(156) 



+ Yld 2 c% -c 2 u z + 2c' 2 cJ<f - 3d 2 c 2 u z + c x c 2 d x q' 



B Matter fluctuations at order O(q ) and 0(P 2 ) 

Unique solution to Eq. (|89|) subject to regularity at the horizon and the boundary 
condition ((73)1 is 

4,o = ln(l - x 2 ) . (157) 

Eqs. (187)) . ()88)1 are second order linear non-homogeneous ordinary differential equa- 
tions which general homogeneous solution can be found analytically Thus it is straight- 
forward to write down a formal solution to ((57)). ((55)1 in quadratures. We do not need 
these explicit expressions for the computation of the hydrodynamic properties of the 
high-temperature cascading gauge theory plasma, so we only argue here that regularity 
at the horizon and the Dirichlet condition at the boundary determine two integration 
constants. We verified numerically that these are the two independent integration 
constants. The latter implies that z^qjZ^ q are uniquely determined. 

Consider first Eq. ((57)) . Its general solution near the horizon x — > + and the 
boundary (1 — x 2 ) = y — ► + takes form 

4o =Ci (2 + In x) + C 2 + C(a; 2 In x) , (158) 



4o =\ (-26 - 44 - jj^j + (4 + 24 + - ^ v + Otf In „) , 

(159) 

where Ci = Ci(Ci,C 2 ) and C 2 = ^(Ci,^)- Regularity at the horizon implies that 
C\ = 0; Dirichlet condition at the boundary further constraints C 2 - 

= 6 + 26 + «fe = 6 (0 - Ca) + 2di (0 ' c *> + • (160) 
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We verified numerically that Eq. ()160)1 indeed has a solution. 

Similarly, the general solution near the horizon x — > + and the boundary (1 — x 2 ) = 
y — > 0+ of Eq. (|HH|) takes form 



-d \nx + C 2 + 0(x 2 \nx) , (161) 



^ 4,o =^ & - ly^ Cx - + y 3 / 2 (4 - 1 & In,) + O (, 2 ) 



(162) 



where we used the power series solution for ip (see Eqs. (|4*Hj) . (JIEJ)) near the horizon 



' ' ^horizon + X 2 ( £ *P horizon + ^ ) + C(x 4 ) , CI 63 ) 



and at the boundary 



^ = -^y + vV /2 -^ 2 + ^ 5/2 ) • ( 164 ) 



Parameters iphorizon, 4> = i>(4>horizon) are uniquely determined from the second boundary 
condition in Eq. 05). In Eq. (gHSD Ci = Cx(Cx,C 2 ) and C 2 = C 2 (Ci,C 2 ). Regularity at 
the horizon implies that C\ = 0; Dirichlet condition at the boundary further constraints 
C 2 : 

= Ci = Ci(0,C 2 ) . (165) 
We verified numerically that Eq. ()165j) indeed has a solution. 



C Matter fluctuations at order ^(q 1 ) and G(P 2 ) 

Unique solution to Eq. (jl(J5J) subject to regularity at the horizon and the boundary 
condition (J74*)l is 

z \ x = ^^(^2x 2 (lna;-2)ln(l-x 2 )+a; 2 dilog(x 2 )+41n(l-x 2 ) N ) . (166) 

As in Appendix El here we only argue that z^z^i solutions to Eqs. (jl()3j) and 
(1104)1 with appropriate boundary conditions exist, and are unique. 

Consider first Eq. ([103)1 . Its general solution near the horizon x —>■ + and the 
boundary (1 — x 2 ) = y — > + takes form 

z 2 ! =d (2 + In x) + C 2 + C(a; 2 In x) , (167) 
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z\ x =i (-2C 2 - 4Ca) + (C 2 + 2di) + y + 0(y 2 hiy) , (168) 

where we used power series solution at the horizon (boundary) for z 2 . Regularity at 
the horizon implies that C\ = 0; Dirichlet condition at the boundary further constraints 
C 2 : 

= C 2 + 2Ci = C 2 (0, C 2 ) + 2Ci (0, C 2 ) • (169) 

We verified numerically that Eq. (J 169)) indeed has a solution. 

Finally, the general solution near the horizon x — > + and the boundary (1 — x 2 ) = 
y -> 0+ of Eq. (fTM)) takes form 

z 2 w>1 =d In x + C 2 + 0(x 2 In x) , (170) 



K* x 4,i =^ Ci - ^ 1/2 Ci + + y^ 2 (c 2 - 1 4 Iny) + O (y 2 ) , (171) 

where we used the power series solutions for ip and z 2 ^ q near the horizon (boundary). 
Regularity at the horizon implies that C\ = 0; Dirichlet condition at the boundary 
further constraints C 2 : 

= Ci = C 1 {0,C 2 ) . (172) 
We verified numerically that Eq. f)172j) indeed has a solution. 
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